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ABSTRACT 


We report on the discovery of FRB 20200120E, a repeating fast radio burst (FRB) 
with low dispersion measure (DM), detected by the Canadian Hydrogen Intensity Map- 
ping Experiment (CHIME)/FRB project. The source DM of 87.82 pc cm ? is the 
lowest recorded from an FRB to date, yet is significantly higher than the maximum 
expected from the Milky Way interstellar medium in this direction (~ 50 pc cm ?). 
We have detected three bursts and one candidate burst from the source over the period 
2020 January-November. The baseband voltage data for the event on 2020 January 20 
enabled a sky localization of the source to within ~ 14 sq. arcmin (90% confidence). 
The FRB localization is close to M81, a spiral galaxy at a distance of 3.6 Mpc. The 
FRB appears on the outskirts of M81 (projected offset ~ 20 kpc) but well inside its 
extended HI and thick disks. We empirically estimate the probability of chance coin- 
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cidence with M81 to be < 107. However, we cannot reject a Milky Way halo origin 
for the FRB. Within the FRB localization region, we find several interesting cataloged 
M81 sources and a radio point source detected in the Very Large Array Sky Survey 
(VLASS). We searched for prompt X-ray counterparts in Swift/BAT and Fermi/ GBM 
data, and for two of the FRB 20200120E bursts, we rule out coincident SGR 1806—20- 
like X-ray bursts. Due to the proximity of FRB 20200120E, future follow-up for prompt 
multi-wavelength counterparts and sub-arcsecond localization could be constraining of 
proposed FRB models. 


Keywords: radio transient sources — fast radio bursts — Milky Way halo 


1. INTRODUCTION 


Fast radio bursts (FRBs) are millisecond-duration, bright radio transients with unknown physical 
origins (Lorimer et al. 2007; Thornton et al. 2013). Although more than 500 bursts have been 
reported thus far (CHIME/FRB collaboration et al., submitted), only 13 published FRBs have been 
sufficiently well localized on the sky to allow host galaxies to be identified. The host galaxies of the 
localized FRBs have redshifts ranging from 0.03 to 0.66, demonstrating the significant distances at 
which FRBs are located!. A small fraction of FRBs have been observed to repeat (Spitler et al. 
2016; CHIME/FRB Collaboration et al. 2019a,b; Fonseca et al. 2020; Kumar et al. 2019), ruling out 
cataclysmic models at least for these sources. Among the repeating FRBs, two have shown periodic 
repetition of as yet unknown origin (CHIME/FRB Collaboration et al. 2020; Rajwade et al. 2020; 
Cruces et al. 2020). 

With this diverse phenomenology, understanding the origins of FRBs is a topical problem in as- 
tronomy. One promising method to unravel FRB physical origins is to study their hosts and local 
environments and compare them with those of the proposed progenitors (Nicholl et al. 2017; Li & 
Zhang 2020). However, FRBs, both repeating and apparently non-repeating, are found in a variety 
of types of host galaxies and it is not yet clear if the two populations are intrinsically different (Bhan- 
dari et al. 2020; Heintz et al. 2020). Therefore, more host localizations are crucial to understand the 
nature of FRB sources. 

The recent discovery of FRB-like radio bursts (The CHIME/FRB Collaboration et al. 2020; Boch- 
enek et al. 2020) and coincident X-ray bursts (Mereghetti et al. 2020; Li et al. 2021; Ridnaia et al. 
2020) from Galactic magnetar SGR, 19354-2154 argue strongly that at least some FRB sources are 
magnetars. Contemporaneous X-ray emission from SGR. 19354-2154 means FRBs may not be solely 
a radio phenomenon. However, to test different proposed FRB progenitor models, multi-wavelength 
follow-up programs are most constraining for nearby sources (« 50 Mpc) due to the sensitivity limita- 
tions of high energy telescopes (e.g. Scholz et al. 2020). Therefore, local universe FRBs are excellent 
candidates to test and constrain different proposed FRB models (for a comprehensive list, see Platts 
et al. 2019, and frbtheorycat.org), and to study in detail the local environments of the FRB sources 
(Marcote et al. 2020; Mannings et al. 2020a; Tendulkar et al. 2020). 

Here we report the discovery of FRB 20200120E, a repeating FRB detected by the CHIME/FRB 
project (CHIME/FRB Collaboration et al. 2018, CHIMEFRBIS hereafter). The FRB dispersion 
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measure (DM) of 87.82 pc cm ? is the smallest reported for an FRB thus far. In 82, we report on 
FRB 20200120E and its three bursts and one candidate burst. In 83, we present evidence against the 
source being associated with the Milky Way or its halo, though the halo origin cannot yet be rejected 
conclusively. We then describe our search for a host and propose that M81, a nearby early-type 
grand-design spiral galaxy at 3.63 + 0.34 Mpc (Freedman et al. 1994; Karachentsev et al. 2002), is 
the most likely FRB host. We further show the FRB is unlikely to be beyond M81. In 84, we discuss 
the implications of this study and conclude in 85. 


2. OBSERVATIONS & ANALYSIS 


(a) (b) 
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Figure 1. Frequency versus time (“waterfall”) plots of three dedispersed bursts detected from FRB 
20200120E. (a) The dedispersed waterfall plots of the three bursts made using total-intensity data, with 
temporal and spectral resolutions of 0.983 ms, and 3.125 MHz, respectively. See Table 1 for the burst prop- 
erties. (b) The waterfall plots of two FRB 20200120E bursts, FRBs 20200120E and 20201129A, made using 
baseband data, which here are binned to have temporal resolution 0.0256 and 0.16384 ms, respectively, and 
spectral resolution 0.391 MHz. The horizontal white lines in the waterfall plots are flagged channels that 
are either bad (due to radio frequency interference) or are missing due to computer nodes in the CHIME 
system being offline. Lastly, the fourth possible burst did not have intensity or baseband data recorded (see 
$2.1), and therefore, not shown here. 


2.1. The Repeating Source FRB 20200120E 


FRB 20200120E was discovered by the Canadian Hydrogen Intensity Mapping Experiment 
(CHIME)/FRB project on 2020 January 20 in the real-time pipeline. The CHIME/FRB instru- 
ment searches 1,024 formed sky beams for dispersed single pulses in the frequency range 400-800 
MHz in total-intensity data, with time resolution 0.983 ms and 16k frequency channels. See CHIME- 
FRB18 for a detailed description of the CHIME/FRB detection system. Offline analysis of the burst 
total-intensity data using fitburst (CHIME/FRB Collaboration et al. 2019a,b,c; Fonseca et al. 
2020) measured a DM of 87.782 + 0.003 pc cm~? (see Tables 1 & 2). 
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This event was also recorded by the CHIME/FRB baseband system, which records full voltage data 
at the Nyquist sampling rate of 800 MHz for all 1,024 CHIME dual polarizarion antennas. Using the 
recorded baseband data and our offline baseband localization pipeline (Michilli et al. 2020), which 
has been calibrated and tested using known positions of Galactic radio pulsars, we have localized the 
FRB to a sky area of z 14 arcmin? (90% confidence region; see Table 2). The baseband localization 
is consistent with that inferred from the multi-beam intensity data detections, but is more precise as 
the full-array baseband data are used to estimate the localization region. 

We detected three bursts and one candidate burst from the FRB source. For the second burst 
on 2020 July 18, the CHIME/FRB baseband system was undergoing maintenance, so we saved 
only total-intensity data at the nominal CHIME/FRB search resolution. On 2020 November 29, 
we detected a third burst for which baseband data were recorded. Following a procedure identical 
to that used for the January 20 burst, we estimate the 90% confidence localization region of the 
burst as: R.A. = 0957™42.81 + 18.89 and Dec = 68?48'57" + 01/39". As the baseband localization 
region of the FRB from the 2020 January 20 burst is completely inside the latter localization region, 
we used the baseband localization region of the FRB in all our follow-up analyses. Later, while 
searching the CHIME/FRB event database, we found a fourth burst detected on 2020 February 6, 
having SNR 10.6, above the 100 threshold that the CHIME/FRB pipeline uses to identify candidate 
FRB events (see CHIMEFRBIS). All the header data metrics, which are discussed in CHIME/FRB 
Collaboration et al. (2018), suggest that the source is astrophysical. However, due to an unidentified 
issue in the FRB data recording system, the intensity data of the 2020 February 6 burst were not 
saved. Therefore, we cannot confirm its astrophysical nature. Hence, we excluded it from the analyses 
in this paper. In Table 1, we report basic properties of this burst from the header data. 

Because the best-fit DMs and sky positions of the three bursts for which we have intensity and/or 
baseband data saved to disk are consistent, we conclude that the FRB is a repeater. Therefore, we 
adopt the Transient Name Server (TNS)? FRB naming convention, which gives the source the name 
of the first detected burst, FRB 20200120E. The TNS names of the bursts on 2020 July 18 and 2020 
November 29 are FRBs 20200718A and 20201129A, respectively. As the intensity data of the FRB 
detected on 2020 February 6 were not saved, we did not request a TNS name for the burst. 


2.2. Burst properties 


Figure 1a shows dedispersed waterfall plots of the three named bursts, made using the estimate 
of each burst’s DM, shown in Table 1, which is optimized to maximize the burst SNR in the total- 
intensity data. We used the calibration methods described by CHIME/FRB Collaboration et al. 
(2019b) and Fonseca et al. (2020) to determine fluences of the three bursts (see Table 1), assuming 
the baseband localization position. Lastly, we employed the same modeling procedures that are 
discussed by CHIME/FRB Collaboration et al. (2019b) and Fonseca et al. (2020) for estimating 
widths, arrival times and scattering timescales from the calibrated total-intensity dynamic spectra 
of the three TNS named bursts. In the total-intensity data, the three named bursts show only one 
component, and therefore, we fitted a single-component profile. 

We also show dedispersed waterfall plots for the two bursts for which we have baseband data saved 
to disk in Figure 1b. We separately optimized DMs of the bursts, FRBs 20200120E and 20201129A, 
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Table 1. Properties of the bursts from FRB 20200120E.* 


TNS Name Day MJD Arrival Time? SNR DM Width Scattering Time Fluence® Peak Flux Density* DME, RM 
(yymmdd) (UTC @ 400 MHz) (pe cm™) (ms) (ms @ 600 MHz) (Jy ms) (Jy) (pe cm^?) (rad m~?) 
20200120E 200120 58868 09:57:35.984(2) 22.9  87.182(3) 0.16(5) «0.23 2.25(12) 1.8(9) 87.789(9) -29.8(5)7 
- 200206° 58885 08:50:45 10.6 88(1) -= — — - — — 


20200718A 200718 59048 22:12:3L882(1) 14.0 87.864 
20201129A/ 201129 59182  13:31:20.8583(6) 19.3 87.812 


5) 0.24(6) «0.17 2.0(7) 1.1(5) - = 


(5 
(4) «01 0.22(3) 2.4(1.4) 1.7(1.2) 87.71(5)  -26.8(3)9 


^ Uncertainties are reported at the 1c confidence level (cl). Reported upper limits are those of the 2c cl. 

b All burst times of arrival are topocentric. 

^ Fluence and peak flux density measurements represent lower bounds as we assumed that the bursts were 
detected at the center of their detection beams. 

4 Optimized DM for the burst detected in the baseband data. 

^ Single beam event with sky position consistent with the FRB baseband localization stated in Table 2. 
However, the intensity data were not saved, so a TNS name has not been assigned. The reported DM and 
timestamp of the FRB were from the header data of the event. 

fBurst parameters were estimated using the total-intensity data. 

9 Both RM measurements have an additional systematic uncertainty of 1.0 rad m7?. 


by aligning sub-structure with the DM. phase? module (Seymour et al. 2019). We estimated best-fit 
DMs of 87.780 + 0.009 and 87.71 + 0.05 pc cm^?, respectively. In the Figure we downsampled the 
data to have temporal and spectral resolution of 0.0256 and 0.16384 ms and 0.391 MHz, respectively. 
The dynamic spectrum of FRB 20200120E reveals downward drifting time-frequency sub-structures 
that are thus far exclusively observed in the dynamic spectra of repeating FRBs (CHIME/FRB 
Collaboration et al. 2019a,b; Hessels et al. 2019; Fonseca et al. 2020; Day et al. 2020). A detailed 
burst analysis of the two bursts using the baseband data is beyond the scope of this paper and will 
be discussed in future work. 


2.3. FRB Rotation Measure 


Following a procedure similar to that outlined by Fonseca et al. (2020), a Faraday rotation measure 
(RM) for FRB 20200120E and FRB 20201129A were measured after applying RM-synthesis (Burn 
1966; Brentjens & de Bruyn 2005) to the burst Stokes Q and U data and detecting a peak in the Fara- 
day dispersion function (FDF). This initial measurement was refined by applying Stokes QU-fitting 
(O'Sullivan et al. 2012), modified to fit simultaneously for parameters characterizing the astrophysi- 
cal signal as well as those corresponding to known systematics. In particular, a delay between the X 
and Y polarizations, arising from different path lengths through the electronics of the system (such as 
cable delay), produces mixing between Stokes U-V parameters (Mckinven et al., in prep.). A best-fit 
model is determined by a nested sampling implementation of QU-fitting that includes an additional 
parameter, 7, characterizing the delay between the two polarizations. Leakage corrected spectra are 
shown in Figure 2 along with model fits. Re-performing RM-synthesis on the cable-delay-corrected 
spectrum results in the FDFs shown in the right panel of Figure 2. We estimate the RM of FRB 
20200120E and FRB 20201129A to be —29.8 + 0.5 + 1.0 rad m^? and —26.8 + 0.3 + 1.0 rad m?. 
Quoted uncertainties correspond to the formal measurement and estimated systematic uncertainties, 
respectively. Ionospheric RM contributions have not been determined here but preliminary analysis 
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Table 2. Major Observables of FRB 20200120E. 


Parameter Value 


R.A.(J2000)^ 09"57™568.7 + 345.6 
Dec. (J2000)* — 68°49/32” + 0124" 


1, b^ 142.°19, +41.°22 

DM*^ 87.818 + 0.007 pc cm? 
DMMw,ne2001 — 40 pe cm? 
DM§iw.ymwis 35 pe cm? 

DMuw, wm 10—40 pc cm~? 

DMIA Na 14—28 pc cm? 

DMi$ rw halo 30 pc cm? 

Max. distance" < 135 Mpc 


RM? —28.3+0.6 + 1.0 rad m? 


^ FRB position determined from baseband data saved for FRB 20200120E using the technique described by 
Michilli et al. (2020). The quoted uncertainty is the 90% cl. 

b Galactic longitude and latitude for the baseband localization central coordinates. 

^ Weighted average DM of the three observed bursts (excluding the 2020 February 6 burst; see Table 1). 

7 Maximum DM model prediction along this line-of-sight for the NE2001 (Cordes & Lazio 2002) and YMW16 
(Yao et al. 2017) Galactic electron density distribution models. 

* DM contribution of the Milky Way assuming that the ISM in the FRB sight-line is dominated by diffuse 
warm ionized medium (WIM); see 83.1. 

f DM contribution of the Milky Way using the Ny — DM relation from He et al. (2013); see 83.1. 

I Milky Way halo prediction from the Dolag et al. (2015) hydrodynamic simulation. The Yamasaki & 
Totani (2020) model predicts a similar value ~ 35 pc cm^?. Note that both these values are smaller than 
the prediction from Prochaska et al. (2019), 50—80 pc cm ?. 

^ Maximum luminosity distance (9096 confidence upper limit) estimated using the Macquart relation (Mac- 
quart et al. 2020); see 83.2. 

i Weighted average RM of the FRBs 20200120E and 20201129A. 


indicates contributions of ~ +(0.2 — 0.4) rad m ? (Sotomayor-Beltran et al. 2013). The difference 
between these two measured RM values therefore is unlikely to be significant. Moreover, both the 
bursts have nearly 10096 linear polarization fraction. The Galactic foreground RM prediction in the 
FRB sight-line from Oppermann et al. (2012) model is —11 +8 rad m ?. The low extragalactic RM 
and DM suggest that the FRB is unlikely to be located in a dense ionized region like a compact H II 
region or star-forming complex (Mitra et al. 2003; Haverkorn 2015; Costa & Spangler 2018; Michilli 
et al. 2018), or in the Galactic center region of a host galaxy (Moss & Shukurov 1996; Krause 2009). 


3. DETERMINING THE DISTANCE AND HOST OF FRB 2020120E 
3.1. Could FRB 20200120E be Galactic? 


In this section, we discuss the possibility that FRB 20200120E is Galactic. The maximum MW 
disk contribution to the DM along the FRB sight-line is 40 pc cm~? from the NE2001 model (Cordes 
& Lazio 2002), or 35 pc cm ? from the YMW16 model (Yao et al. 2017). The observed FRB DM is 
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Figure 2. Summary plots of the RM detection methods of QU-fitting and RM-synthesis applied to the 
cable delay corrected spectra of FRB 20200120E (top row) and FRB 20201129A (bottom row). Left panel: 
Stokes Q, U normalised by the total linear polarization (P = y Q? + U?) and polarization angle, v, as a 
function of frequency with corresponding model fits. Frequency channels with significant polarized signal 
are highlighted through a greyscale that saturates at higher S/N. Right panel: The results of RM-synthesis 
showing each event’s FDF constrained near the peak. 


significantly larger than the DM predictions of either model even after taking into account an 2096 
systematic uncertainty (Cordes & Lazio 2002; Yao et al. 2017). If the FRB is Galactic, an H I! region 
could contribute to the DM-excess of the FRB as discussed by Patel et al. (2018). We checked the 
Anderson et al. (2014) Hit region catalog, which is claimed to be complete for Galactic H II regions 
other than large diffuse and young hypercompact H II regions, and found none, unsurprising given 
the high Galactic latitude of the source (b = +41.°22; see Table 2). We also did not find any CO 
emission within the FRB localization region in the Planck all-sky CO map (Planck Collaboration 
et al. 2014) ruling out the presence of a young hypercompact H 1 region in a dense CO clumps of 
a molecular cloud complex (Dame et al. 2001). Finally, we searched the FRB field-of-view in the 
Wisconsin Ha Mapper Northern Sky Survey (WHAM; Haffner et al. 2003) for the presence of any 
extended and diffuse Ha excess clump but did not find any. Therefore, it seems unlikely that an H 
II region is responsible for the observed DM-excess. 

We consider two additional, independent maximum Galactic DM estimates. Using the observed 
Galactic hydrogen column density Ny in the FRB sight-line from the HI 47 survey (HI4PI) (HIAPI 
Collaboration et al. 2016), we find Ng = 5.8 x 10? cm~?. Using an empirically derived Ny — DM 
relation from He et al. (2013), we calculated the Milky Way contribution to the DM in the direction 
of the FRB to be ~14—28 pc cm ^, significantly smaller than the DM of the FRB. Note, however, 
that this Ny—DM relation was estimated using radio pulsars, which are generally located near to the 
Galactic plane. 
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We also independently estimated the Milky Way disk DM that, at high Galactic latitudes, is 
dominated by the warm ionized medium (WIM) extending to a vertical distance ~ 1—2 kpc above 
the Galactic plane (e.g., Hill et al. 2015). From the Planck all-sky free-free emission map (Adam 
et al. 2016), we find a total emission measure (EM) in the FRB sight-line of ~ 7.8 pc cm ?. To 
estimate the WIM-dominated MW disk DM, we used a range of free electron density (Ne) values of 
the WIM from the literature (Reynolds et al. 1995; Reynolds 2004; Gaensler et al. 2008; Velusamy 
et al. 2012; Ocker et al. 2020), Me zz 0.2—0.9 cm ?, and computed DM = EM/N, ~ 10—40 pc cm ?. 
Moreover, Hill et al. (2007) estimated the mean dispersion measure, DM sin |b|, for high-EM sight- 
lines (EM sin |b| > 2 pc cm 9) through the WIM to be 14.8 0.9 pc cm ?. At this FRB’s Galactic 
latitude, this relation gives a mean Galactic DM of 22.5 + 1.4 pc em ?. From all these estimates, it 
seems highly unlikely the source is within the Milky Way disk. 

The DM contribution of the MW halo, DMmw halo, which consists primarily of hot ionized circum- 
galactic gas extending out to a Galactocentric radius of ~ 200 kpc, is poorly constrained (Keating & 
Pen 2020). If we assume a halo DM contribution of 50—80 pc cm ? as proposed by Prochaska et al. 
(2019), the FRB could be within the MW halo. However, two other MW halo DM models, those of 
Yamasaki & Totani (2020) and Dolag et al. (2015), predict DMMW halo ~ 30 pc cm ?, which support 
an extragalactic origin for the FRB. 

If the FRB is a MW halo object, the DM-excess of ~ 50 pc cm"? must be contributed by the 
hot coronal gas in the MW halo. Using different tracers of the hot gaseous medium (10? — 10’ K; 
'Tumlinson et al. 2017) that include emission and absorption lines of highly ionized species in the UV 
and X-ray (Miller & Bregman 2015) and constraints from the diffuse soft X-ray background (Henley 
& Shelton 2013), the electron density of the MW halo is estimated to be ~ 10^? cm ? close to the 
Milky Way disk (Bregman & Lloyd-Davies 2007) and ~ 1074 cm ? at ~ 50 kpc (Bregman et al. 
2018), further reducing to ~ 107° cm ? near the virial radius of the Milky Way (Kaaret et al. 2020). 
Assuming a gas filling factor of unity and electron density of the coronal gas, Ne ~ 1078—1074 cm^?, 
we estimate the distance to the FRB source to be ^ 50—500 kpc. As the Milky Way virial radius is 
~ 200 kpe (Dehnen et al. 2006), finding the FRB source at a distance > 200 kpc implies that the 
source is extragalactic. Therefore, in further discussion, we use 50-200 kpc as a plausible distance 
range for a halo object. 

In this scenario, it is possible that the FRB source could be associated with either a Milky Way 
satellite galaxy (Karachentsev & Kaisina 2019; Kaisina et al. 2019) or globular cluster (Harris 2010; 
Gaia Collaboration et al. 2018; Vasiliev 2019), but no such cataloged source exists within or near 
(within ~ 10° radius circular sky area) to the FRB localization region (Contenta et al. 2017; Simon 
2019). 

If the source is in the MW halo, the observed bursts could correspond to super-giant pulses (SGPs) 
from one of a young neutron star, millisecond pulsar, or rotating radio transient (RRAT). We now 
consider each of these possibilities in turn. 

It is difficult to explain the existence of a young neutron star in the halo, given that these objects 
are expected to form only near the Galactic Plane. The fastest known runaway OB stars and young 
pulsars both have space velocities ~ 1000 km s^! (Brown 2015; Du et al. 2018; Chatterjee et al. 
2005), which would require ~ 108 yr to traverse 50 — 200 kpc. This is significantly longer than the 
typical lifetime of an OB star (tens of Myrs; Crowther 2012) and inconsistent with the expectation 
of a young neutron star. On the other hand, it has been suggested (e.g. Giacomazzo & Perna 2013) 


9 


that young neutron stars can be formed by compact object mergers, which could occur in the halo. 
However, this formation mechanism has not been confirmed as actually occurring in nature. 

Few millisecond pulsars (MSPs) are also known to emit SGPs (Johnston & Romani 2004). In 
principle, an isolated MSP can exist at a distance of ~ 50 kpc. However, Galactic isolated MSPs 
that are known to produce SGPs are rare; only two (PSRs B1937+21 and B1821-24: Cognard et al. 
1996; Romani & Johnston 2001) such sources have been seen to do so out of a sample of ~ 450 
known MSPs (Manchester et al. 2005). Moreover, assuming a SGP duration of ~ ps, the maximum 
isotropic energy observed from the brightest SGPs of PSR B1937+21 is ~ 107° erg (McKee et al. 
2019). The isotropic burst energy of the FRB source, if located at a distance of ~ 50 kpc, would 
be ~ 10”? erg, 100 times brighter than the brightest giant pulses observed from PSR B1937+21. 
Similar analysis using PSR B1821—24 SGPs would give an even larger energy difference (Johnson 
et al. 2013). Note that the energy difference would further increase if we assume the FRB source 
distance > 50 kpc. 

Alternatively, the FRB source could potentially be explained as Crab-like SGPs from a ~ 10? year 
old rotating radio transient (RRAT) in the halo. However, it is unclear if an old RRAT can produce 
such pulses. For instance, the RRAT with the largest isotropic energy estimated using the distance 
and mean flux density values from the RRAT catalogue* is RRAT J1819—1458 with an energy of 
~ 10” erg. This is again two orders of magnitude below the energy needed to power FRB 20200120E 
at a distance of 50 kpc. 

Finally, FRB 20200120E shows complex spectral and temporal downward-drifting sub-structures 
(see Figure 1b), which have previously been established as a characteristic spectro-temporal feature 
of repeating FRBs (CHIME/FRB Collaboration et al. 2019b; Hessels et al. 2019; Fonseca et al. 2020; 
Day et al. 2020). Such structures are not known to be seen in pulsar or RRAT spectra, although 
some Crab SGPs have shown similar complex structure (Hankins & Eilek 2007). 

Overall, if FRB 20200120E is a MW halo object, we conclude that it would be the most distant 
Galactic neutron star yet discovered and would also need to be unusually energetic compared to 
known objects. 


3.2. Host galaxy search 


In light of the challenges faced when associating FRB 20200120E with an object in the MW halo, 
we next consider whether it could be associated with an external galaxy. If we assume DMhaio = 30 pc 
cm ^? as predicted by the Dolag et al. (2015) and Yamasaki & Totani (2020) models, the extragalactic 
DM of the FRB, DMgg, is 18 and 23 pc cm ? for the NE2001 and YMY16 Milky Way DM models, 
respectively (see Table 2). For a negligible host DM contribution, we estimate the maximum redshift 
of the FRB to be Zmax 7 0.03 (90% confidence upper limit), or maximum luminosity distance of 
135 Mpc, using the Macquart DM—z relation (Macquart et al. 2020). Therefore, if the FRB is 
extragalactic, we expect a nearby host galaxy within its location error region. 

One possibility is that there is a faint dwarf galaxy, like that of FRB 121102, in the localization 
region of FRB 20200120E. An FRB 121102-like star-forming dwarf galaxy (M, — —17 AB mag; 
Tendulkar et al. 2017), if located at Zmax, would have r-band magnitude ~ 19 AB mag. Within the 
FRB 9096 confidence localization region, shown as a red ellipse in Figure 3, we identify one galaxy 
in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Survey photometric catalogue 


4 http://astro.phys.wvu.edu/rratalog/ (visited on 19/12/2020) 
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Figure 3. The Digital Sky Survey (DSS) RGB image of the region around M81. The red ellipse represents 
the 90% confidence localization region of FRB 20200120E. Source 1 is the cataloged M81 H 1 region, 
[PWK2012] 31 (Patterson et al. 2012), Source 2 is an X-ray source, [SPZ2011] 8 (Sell et al. 2011), Source 3 
is an M81 globular cluster, [PR95] 30244, and Source 4 is the VLASS point radio source, VLASSIQLCIR 
J095756.104-684833.3 (Gordon et al. 2020). All the sources are found in the outer disk of M81. The inset 
image is the 21-cm line view of the M81 circumgalactic medium (Chynoweth et al. 2008); the dashed magenta 
box is the DSS image field-of-view. 


(Dey et al. 2019), 2MASX J09575586--6848551 with r-band magnitude — 15.35 AB mag. However, 


it is known to be at z = 0.19395(2) which is significantly further away than Zmax (Huchra et al. 2012). 


11 


Therefore, any dwarf host galaxy within the 90% confidence region must be fainter than the FRB 
121102 host. 

More interestingly, we find the FRB sight-line has a sky-offset from the M81 center, a nearby 
grand-design spiral galaxy, of 19.6. At the 3.6 Mpc distance of M81 (Karachentsev et al. 2002), this 
sky offset corresponds to a projected distance of 2073 kpc? from the center of M81, well within the 
extended Hi disk (see Figure 3) and thick disk of M81 (~ 25 kpc; Tikhonov et al. 2005). The FRB 
localization region is in the location where tidal interaction among M81 group members has resulted 
in the formation of star forming clumps. Observations in the visible band give only marginal sign 
of these sources, but they are extensively studied in the radio and X-ray bands. These sources are 
discussed in 83.5. Moreover, between the M81 group and Milky Way, we did not find any cataloged 
field or Milky Way halo satellite galaxy. All these observations make M81 a plausible host galaxy for 
FRB 20200120E. 


3.3. Can the proximity to M81 be by chance? 


We now estimate the chance coincidence probability (Pec) of finding an M81-like bright galaxy close 
to the FRB localization region. We define Pec = Aca/Acume, where Aga is the total angular sky 
area spanned by MSI-like or brighter galaxies that are visible to CHIME, and Acyrme is the total sky 
area visible to CHIME (sky area above Dec = —10°; 61% of the total sky area). To be conservative, 
we remove the Milky Way sightlines where the DM-excess of FRB 20200120E is less than the ~ 10% 
systematic error on the maximum of the two different Milky Way DM model predictions (Cordes 
& Lazio 2002; Yao et al. 2017), which we define as DMex. This DM-excess constraint (hereafter 
C1) removes 10% of the total sky visible to CHIME (mostly consisting of the Galactic plane), and 
we estimate Acyrmm = 20600 deg?. Note that the CHIME sensitivity changes with declination, but 
this effect is likely insignificant in our case as all the nearby bright galaxies are within 30? of the 
zenith in the CHIME primary beam. Next, we use the catalog of the local volume galaxies?, which 
is complete for M81-like bright galaxies, and find three galaxies other than M81 that have extinction 
corrected B-band magnitudes, mg < 7.5, mg of M81': M31 (Andromeda, mg = 3.7 at 770 kpc), 
M33 (Triangulum, mg = 6.1 at 930 kpc), and IC 342 (mg = 7.2 at 3.28 Mpc). To estimate the 
total sky area of these galaxies, we use a circular region with an angular radius equivalent to a 20 
kpc projected offset (of FRB 20200120E from M81) at their respective distances, which are e 1.°49, 
1.°23, 0.°35, and 0.°31 for M31, M33, IC 342, and M81, respectively. Using these values we estimate 
Acai 7: 12.4 deg, and hence, Pee = 6 x 1074. 

As the presence of M81 is inferred post-hoc, it is essential to correct the chance coincidence prob- 
ability for the problem of multiple testing (also known as the look-elsewhere effect), that tends to 
increase the false-positive rate of a discovery (Type I error; Maxwell et al. 2017). To account for this, 
we use the Bonferroni correction procedure (Bender & Lange 2001; Armstrong 2014). We consider all 
FRBs in the first CHIME/FRB catalog (CHIME/FRB Collaboration et al., submitted) that satisfy 
the excess-DM constraint C1. The Bonferroni correction inflates the Pee to 1 — (1—6 x 10~+)NFRB.DMex , 
where Nprp.pm., is the number of FRBs in the CHIME/FRB catalog with the DM-excess € DM. 
Figure 3.3 shows the Pece as a function of DM,,. As FRB 20200120E has the lowest DM,, in our 
sample, Nprp7s = 1, and therefore, Pee = 6 x10 *. To be conservative, we also count the num- 


? 9096 c.l 
8 https:/ /www.sao.ru/lv/lvgdb/introduction.php (visited on 19/12/2020) 


7 Though we have used B band magnitudes that are provided by the catalog, the results would not change if we use 


other optical band magnitudes. 
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captionThe chance coincidence probability of finding an M81-like galaxy as a function of the DM-excess 
(DM) of FRBs in the first CHIME/FRB catalog (submitted) that satisfy the DM-excess constraint, C1, 
as discussed in 83.3. This takes into account the look-elsewhere effect by incorporating the number of 
FRBs in the CHIME/FRB catalog with DM-excess < DMex as a trial factor. FRB 20200120E has the 
lowest DM, which gives Pee = 6 x 1074. At the DM of FRB 20200120E = 87.82 pc cm^?, Po. = 0.7%, 
which we consider as our conservative Pec estimate. Moreover, the results would not change if we use 
minimum of the two different Milky Way DM model predictions (Cordes & Lazio 2002; Yao et al. 2017) to 
estimate DM, (instead of the maximum predicted values used in the DM-excess constraint C1). Lastly, 
we also showed the chance coincidence probability when only Catalog-1 FRBs with saved baseband data 
are used in correcting for the multiple testing problem. 


ber of CHIME FRBs with DM-excess < 87.82 pc cm^?, the DM of FRB 20200120E, and estimate 
Nerrp,s7.s2 = 11. This gives a value of Pee = 0.007. 

There are several factors that make our Pec estimate conservative. First, the Bonferroni correction 
becomes overly conservative as the number of events increase, and hence, undermines the significance 
of an unlikely observation (Perneger 1998; Nakagawa 2004). Second, we search prospective hosts 
for only those FRBs that have saved baseband data, which constitute a small fraction of the first 
CHIME/FRB catalog FRBs. Therefore, it could be argued that we should have used only them 
when accounting for the multiple testing problem (Streiner & Norman 2011; Maxwell et al. 2017). 
We also considered other nearby galaxies, except M81 satellite galaxies which are discussed in §3.5, 
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that have projected angular offset less than or equal to that of M81 (19.'6) from the FRB, and for 
none is Pee < 10%, even before correcting for the multiple testing problem. However, note that the 
formalism of chance coincidence probability favors brighter galaxies over the fainter ones. If FRBs 
preferentially originate in a specific galaxy type, then the Pec is not a good proxy of true a host 
association probability. We need more host associations to test the latter possibility. 

Bayesian hypothesis testing is another method proposed to avoid the problem of multiple testing 
(Scott & Berger 2006; Gelman et al. 2012). However, its success in avoiding the problem of mul- 
tiple testing strongly depends on the choice of priors (Rouder 2014; de Heide & Grünwald 2020). 
With better knowledge of the nature of FRB hosts, it will eventually be possible to use a Bayesian 
framework to estimate true host association probabilities. 


3.4. Could FRB 20200120E lie beyond M81? 


Here we discuss the possibility of the FRB source being located beyond M81. In such a scenario, 
the CGM of M81 and its neighboring galaxies will also contribute to the FRB DM. 

M81 is a part of the nearby “M81 group,” which contains prominent galaxies like M82, NGC 2403, 
NGC 4236, and several dwarfs (Karachentsev et al. 2002). Some of its members are in the process 
of merging, which makes the M81 group CGM rich in metals and gas (Al Najm et al. 2016). We 
estimate the DM contribution of the M81 halo assuming that the true host lies beyond M81. 

The FRB 20200120E sight-line passes through the M81 halo with a very small impact parameter 
(~ 20 kpc), which is much less than the M81 virial radius of ~ 210 kpc (Oehm et al. 2017). Therefore, 
we expect the M81 halo to contribute considerably to the FRB DM if the FRB is located beyond 
M81. To quantify this effect, we assume a density profile and several major M81 halo parameters, 
such as halo mass and virial radius. We first consider two halo gas profiles, proposed by Maller & 
Bullock (2004) (hereafter MB04) and Mathews & Prochaska (2017) (hereafter MP17)*. Both profiles 
are modified versions of the Navarro-Frenk-White (NFW) density profile (Merritt et al. 2005). For 
both profiles, we assume the following parametric form (Prochaska et al. 2019): 


Po 

bot E 
Here for the MB04 profile, we adopt a = 2, yo = 4, and for the MP17 profile, a = 2, yo = 2. In 
Equation (1), po is the central halo density, and Rs = Roo9/c, where Roo is the scale radius that 
encloses a density of 200 times the critical density of the Universe, pai, = 3H?/87G, and c is the 
concentration parameter defined as the ratio between the virial and scale radius of a halo. T'he values 
of these parameters for the M81 halo are c = 10.29, Ro = 210 kpc, and pọ = 7.21 x 10? Ma pc, 
taken from Oehm et al. (2017). As suggested by Prochaska et al. (2019), we terminate the density 
profile at the M81 virial radius. Additionally, we need to assume a value for the fraction of baryons 
that is retained inside the virial radius of the M81 halo, fy halo: For this, we consider two values: (1) 
fy.naio = 0.75 that assumes ~ 25% of the baryons exist in the galaxy as the interstellar medium (ISM), 
stars, and compact remnants (Fukugita et al. 1998), and (2) fp,halo = 0.40 which is a lower limit that 
Hafen et al. (2019) found for  10!?M, mass halos in the Feedback in Realistic Environment (FIRE) 
simulation. The FIRE simulation is well suited to study the CGM of simulated galaxies with high 


8 The Dolag et al. (2015) and Yamasaki & Totani (2020) models are tailored to the Milky Way, and therefore, cannot 


be used here. 
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Figure 4. M81 DMhajo as a function of impact parameter of the FRB assuming M81 is a foreground galaxy. 
We used two modified versions of the Navarro-Frenk-White (NFW) halo profile discussed by Prochaska 
et al. (2019): MB04 (a = 2, yo = 4) and MP17 (a = 2, yo = 2); see 83.4 for detailed description. The 
DM profiles are plotted for two M81 halo baryon fractions: (1) fy nai = 0.75 (dash lines), a fiducial value 
that presumes that the halo has retained the mean cosmic baryons and that 2596 of these baryons are in 
the galaxy as stars, collapsed objects, and interstellar medium (ISM; Fukugita et al. 1998);(2) fb halo = 0.4 
(solid lines), the minimum value that Hafen et al. (2019) found in the FIRE simulation for a halo of mass 
^ 101? Mo. The region between the vertical dotted lines represents the range of the FRB projected distance 
from the center of M81 given the uncertainties in the FRB 90% confidence localization region. As can be 
seen from the plot, even the most conservative scenario (MB04 profile with fp halo = 0.4; solid purple line) 
predicts DMpaio greater than the DM-excess of FRB 20200120E, 18 and 23 pc cm-?, using NE2001 and 
YMW16 model predictions, respectively. This suggests that it is unlikely that the FRB host is beyond M81. 
We also consider a non-NFW halo profile proposed by Pen (1999) (P99) with two core radius (re) values 
from Keating & Pen (2020): re = R200 and re = 0.86 R200, where Roo is the virial radius of the M81 halo; 
the large value of re allows more halo gas to be expelled from the M81 halo. Even in this scenario, we find 
the M81 DMhaio predicted by P99 model is significantly larger than the DM-excess of the FRB (see 83.2). 
Here we have included only the M81 halo contribution; adding contributions from other members of the 
M81 group would strengthen our conclusion. Note that if the FRB is associated with M81, the M81 DMhajo 
contribution would be likely smaller than the estimated values (< 50%). 
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resolution and also takes into account the effect of gas inflow and outflow along with other factors 
that play important roles in the evolution of galaxies (Hopkins et al. 2014, 2018). To convert the 
dark matter density to that of baryons, py, we use the cosmic baryon fraction = 0Q,/Qm = 0.158 
(Ade et al. 2016). Finally, we estimate the free-electron density using the relation: n, = f.(pp/mp)’, 
where fe is the number ratio between free electrons and baryons in the halo (~ 7/8; Tumlinson et al. 
2017), and m, is the proton mass. Using these parameters, we estimate the M81 DMhpaio using the 
following equation: 


roux Ri 
DMhalo = 2 f nedl. (2) 
0 


In Equation (2), rua is the maximum radius of integration through the M81 halo (Rgoo in our 
analysis) and R, is the FRB impact parameter. 

Apart from the modified NFW-halo density profiles, we also consider the entropy-floor singular 
isothermal sphere model of Pen (1999). The model invokes two phases of halo gas where in the inner 
region the gas is heated to constant entropy, and at radius (r) > core radius (re), the gas traces the 
halo mass isothermally. The gas density profile is defined as : 


2 ls 
pr) = (42) jl -- E Es 1.59. (3) 
T zl + gsln( s)? &if r < re 
Here f, = 0.06h71° and Veire = (10GHoMp)! ~ 154 km s^! are the gas fraction and circular velocity 
of the M81 halo of mass Mp = 1.3 x 10? Mo (Oehm et al. 2017) taken from Pen (1999) and Mo et al. 
(1998), respectively. We use two values of r, in our analysis as suggested by Keating & Pen (2020): 
re = Rago and re = 0.86 Rao. To estimate the M81 halo DM, we use the procedure as discussed 
above. 

Figure 4 shows our estimates of the M81 DM),),. as a function of R,. At the FRB R, ~ 20 kpc, 
the minimum DM halo value (estimated using the MB04 profile and fb halo = 0.4) is larger than the 
FRB DM-excess of 18—23 pc cm ?, see Table 2). This argues against the FRB being beyond or even 
on the far side of the halo. To check if our most constraining M81 DMmyg, estimate is conservative, 
we estimated the Milky Way DMmgay using the MB04 profile, fb haio = 0.4, and other Milky Way 
halo parameters from Prochaska & Zheng (2019), and found DMgys z 13 pe cm ?. This is smaller 
than the Dolag et al. (2015) and Yamasaki & Totani (2020) predictions in the FRB sight-line, ~ 30 
pe cm~’. This demonstrates that our choice of halo density model does not bias the M81 DMhato 
analysis. 

Note that in our analysis, we consider only the M81 halo contribution. However, other members of 
the M81 group would also contribute to the FRB DM. For instance, we can calculate DM), of the 
M81 group by only considering the two nearest massive M81 group members other than M81, NGC 
3077 (Rx ~ 36 kpc) and M82 (R, ~ 59 kpc), with their respective halo parameters from Oehm 
et al. (2017). The total DMnaio of the M81 group we estimate using the MB04 density profile with fi 
= 0.4, the most conservative case in our analysis, is ~ 75 pc cm ? > DMgg = 53—48 pc cm ? (for 
a negligible Milky Way halo contribution). This further strengthens the conclusion that the FRB is 


9? We presume a flat ACDM (Planck Collaboration et al. 2016) model with the matter density, 
Om = 0.308, baryonic matter density, €, = 0.0486, dark energy density, Qa = 0.691, and Hubble constant, 
Ho = 100 h km s-! Mpc^! with h = 0.6774. 
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unlikely to be beyond the M81 group. Should the FRB source turn out to be behind M81, it would 
suggest that the M81 group has lost most of its halo baryons. 

As shown in Figure 3, the FRB localization region contains a significant amount of HI gas. Tikhonov 
et al. (2005) argued that the M81 thick disk extends to a galactocentric radius of 25 kpc, in which 
case the observed HI flux should be considered as a part of the M81 interstellar medium and hence, 
should make an additional DM contribution. In order to estimate the contribution of the Hı disk 
(DMg), we use the Ny—DM relation from He et al. (2013). As discussed in 83.1, this relation is 
derived using nearby Milky Way radio pulsars, and so it is unclear whether it is valid for the M81 HI 
disk. Noting the similarity between the Milky Way and M81 projected Hı distribution (Westpfahl 
et al. 1999), we use this relation for a rough estimate of DMy;. In the Very Large Array (VLA) 
image of the M81 group created by de Blok et al. (2018), the mean integrated HI flux intensity 
within the 90% localization region is 0.21 Jy beam! x km s^!. The VLA D-configuration beam 
is modelled as an ellipse with major (0,) and minor (6) axes 38."01 and 30."91, respectively, and 
position angle — 75?.5 that we get from the header data of the 21-cm map fits file made by de Blok 
et al. (2018)!°. Assuming the HI gas to be optically thin, we estimate the integrated hydrogen column 
density, Ny & 2 x 107° cm?. 

Using this Ny value in the Ny—DM relation, we estimate DMy, ~ 5 — 10 pc em^?, over and above 
the contribution from the halo, further diminishing the probability that the FRB lies beyond the 
M81 group. 

We conclude that FRB 20200120E is unlikely to lie beyond the M81 group. Therefore, if the FRB 
is an extragalactic source, its host is most likely located in the M81 group. 


3.5. Interesting M61 Group Sources 


In 83.4, we have shown that the FRB source is unlikely to be located beyond M81. Moreover, if 
we consider the MP17 halo density profile, the FRB is most likely to be located at the near side or 
in front of the M81 halo. However, with the MB04 and P99 halo density profiles, it is possible for 
the FRB to be present in the M81 extended Hı disk. In any case, we searched for any cataloged 
M81 group satellite galaxy and did not find any within the FRB localization region. Within the M81 
halo of sky-radius ~ 3.°3 (projected angular offset corresponding to the M81 virial radius, 210 kpc; 
Oehm et al. 2017), we found 42 dwarf satellite galaxies from the literature (Van den Bergh 1966; 
Karachentseva et al. 1985; Karachentseva 1968; Boerngen & Karachentseva 1982; Karachentseva 
et al. 1985; Caldwell et al. 1998; Froebrich & Meusinger 2000; Boyce et al. 2001; Karachentsev & 
Karachentseva 2004; Chiboucas et al. 2009; Smercina et al. 2017; Okamoto et al. 2019). We found 
a young dwarf irregular galaxy, Holmberg IX, and a dwarf spheroidal galaxy, KDG 64, at offsets 
of 13/9 and 14/3, respectively, from the center of the FRB localization region. Note that the two 
galaxies have the offset greater than the 10 times their half light radius. However, using the M81 
satellite number density, pg ~ 1.2 deg ?, and assuming a Poisson distribution within the M81 halo 
as found in the semi-analytic and N-body gas dynamics studies by Kravtsov et al. (2004), the chance 
coincidence probability (Pec) of finding these two galaxies near the FRB localization region can be 
estimated using the relation from Bloom et al. (2002), Pec = 1 - es^, where A is the circular sky 
region of radius equals to the angular offset of the galaxies. We find Pee > 10% for the galaxies. 


10 The image data can be directly downloaded from https://www.astron.nl/~blok/M81data/. 
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Within or close to the FRB localization region, we found an M81 globular cluster, [PR95] 30244 
(Source 3 in Figure 3; Nantais & Huchra 2010). Though unremarkable, its presence in the FRB 
localization region is still noteworthy as globular clusters in the Milky Way are known to host several 
exotic systems, like millisecond pulsars, blue stragglers and low-mass X-ray binaries (LMXBs) (Brodie 
& Strader 2006). Moreover, as shown in Figure 3, we found an H 11 region, [PWK2012] 31 (Source 1), 
and an X-ray source, [SPZ2011] 8 (Source 2), within the FRB localization region. The presence of an 
HII region confirms that in situ star formation is actively taking place within the extended H1 M81 
disk. Additionally, Mouhcine & Ibata (2009) found a number of young stellar systems with ages of 
only a few tens of Myr in the M81 extended Hı disk, which were formed during the tidal interaction 
of M81 with surrounding companion galaxies, most prominently NGC 3077 and M82. Therefore, 
it seems possible that a young neutron star — a possible FRB counterpart — associated with M81 
is present within the FRB localization region. Lastly, Sell et al. (2011) estimated the counts in 
different X-ray sub-bands in the range 0.5—8 KeV. Using the X-ray color-based source classification 
first proposed by Prestwich et al. (2003), if the X-ray source [SPZ2011] 8 is indeed an M81 source, we 
find that it can be either a high mass X-ray binary (HMXB), LMXB, or thermal supernova remnant. 
A future more precise FRB localization will tell us if the FRB is actually associated with any of these 
M81 sources. 


3.6. Search for a persistent radio source 


We also searched for a persistent radio source, like the one seen coincident with FRB 121102 by 
Chatterjee et al. (2017), within the 90% confidence localization region of the FRB. We searched 
archival data of the following surveys: the NRAO VLA Sky Survey (NVSS; Condon et al. 1998), 
the VLA Sky Survey (VLASS; Lacy et al. 2016), the Westerbork Northern Sky Survey (WENSS; 
hengelink et al. 1997), and the Tata Institute of Fundamental Research Giant Metrewave Ra- 
dio Telescope Sky Survey (TGSS; Intema et al. 2017). From the search, we identified only one 
point source in the VLASS image (2—4 GHz frequency band) above the 5e noise threshold (~ 0.6 
mJy beam 1). In the Canadian Initiative for Radio Astronomy Data Analysis VLASS Epoch 1 Quick 
Look (CIRADA VLASS QL) Catalog (Gordon et al. 2020), this source is cataloged as VLASSIQLCIR 
J095756.10+684833.3 with R.A = 09^57"56* and Dec.— +68748™33° (Source 4 in Figure 3). The 
estimated integrated and peak flux density of the source are 1.27 + 0.19 mJy and 1.35 + 0.11 
mJy beam", respectively. Using the 5-o upper-limit from the NVSS image (1.4 GHz) of the FRB 
field-of-view of 2.15 mJy beam"! and assuming a power-law dependence of the source flux density 
i.e., S, x S°, we estimated a lower limit a > —0.61. Note that the VLASS calibrated images are 
known to have a systematic that underestimates the flux density of identified sources by ~ 10% 
(Gordon et al. 2020). Incorporating this in our spectral index estimate would increase the derived 
lower limit. Active galactic nuclei (AGN) are known to show variable light curves at all frequencies. 
Given a ~% 20-yr gap between the NVSS and VLASS observations, the derived spectral index upper 
limit is not constraining if the radio source turns out to be a background AGN. 

We then compared the VLASS source with the one that was found spatially associated with FRB 
121102 (Chatterjee et al. 2017). At 3.6 Mpc, the isotropic luminosity of the radio source at 3 GHz, 
~ 10? erg s !, would be ~ 10^ times smaller than that of the FRB 121102 radio source. The 
isotropic luminosity would be ~ 105 — 10’ times smaller if the radio source is at distance 50—200 
kpc. We can rule out a canonical stellar-mass X-ray binary, as such a source always has radio 
luminosity smaller than 1033 erg s! at ~ 1 GHz even when flaring (Gallo et al. 2018; Reines et al. 
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2020). A pulsar wind nebula (PWN) can give rise to a flat spectrum radio source (a > —0.4; Reynolds 
et al. 2017; Slane 2017). We searched for a possible optical counterpart of the radio source in the 
CFHT/MegaCAM image from Chiboucas et al. (2009) and found none. The 5-o of r-band limit = 
25.65 AB mag corresponds to an r-band flux of 1.5 x 1099 erg s^! at 3.6 Mpc. If the radio source is a 
Crab-like PWN with radio-to-optical luminosity ratio < 10^? (Volpi et al. 2008), it should have been 
detected in the CFHT/MegaCAM r-band image. Moreover, Sell et al. (2011) observed M81, including 
the FRB field-of-view, using the Chandra X-ray Observatory and identified 276 X-ray point sources 
with sensitivity ~ 10°” erg s*!. However, there is no X-ray point source spatially coincident with the 
radio source. If the radio source is an M81 supernova remnant or PWN, its X-ray/radio flux ratio 
is likely < 10~?; otherwise, it is most likely a background AGN. Note that these constraints are also 
valid if the FRB source is located at 50—200 kpc. The persistent radio emission from X-ray binaries 
are orders of magnitude lower than their X-ray luminosities (Koljonen & Russell 2019), therefore, 
a Galactic X-ray binary as the counterpart of the persistent radio source is also disfavoured from 
the constraints discussed above. Future follow-up observations of the persistent radio source and a 
milliarcsecond localization of the FRB may one day tell us if the radio source is associated with the 
FRB. 


4. DISCUSSION 


We have shown that the sky location of the low-DM repeating FRB 20200120E appears super- 
imposed on the extended HI and thick disks of the nearby spiral galaxy M81. Moreover, the low 
DM-excess of the FRB suggests that the FRB source is unlikely to be located beyond M81 group 
(~4 Mpc). We searched for galaxies closer than those associated with the M81 group within or 
near to the FRB localization region in the catalog of the local volume galaxies (Karachentsev et al. 
2013) and found none. Additionally, the coincidence probability of finding an M81-like galaxy close 
to the FRB localization region is small (« 1%). Therefore, if extragalactic, the FRB is most likely 
associated with M81 which would make it by far the closest extragalactic FRB yet known. Lastly, 
given the observational constraints, we cannot reject the Galactic origin of the FRB. 


4.1. Constraints on the Milky Way halo DM contribution 


Under the assumption that FRB 20200120E is extragalactic, it can be used to set an upper limit 
on the MW halo DM contribution in this direction. If we consider the lowest of the two MW DM 
model estimates, DMmw,ymw16 = 35 pc cm? and conservatively assume negligible IGM and host DM 
contribution, then DMyw halo < 53 pc cm~’. This would be inconsistent with most of the DMmw halo 
phase space proposed by Prochaska et al. (2019), i.e., DMyai; = 50—80 pc cm ?. On the other hand, 
both the Dolag et al. (2015) and Yamasaki & Totani (2020) models predict DMy4 ~ 30 pc em ?, 
lower than our upper limits. The halo may be clumpy (Kaaret et al. 2020; Keating & Pen 2020), so 
it may still be possible to have significant variations in DMpa). along different sight-lines. More such 
low-DM FRB localizations will help in constraining the structure and composition of the MW halo. 


4.2. Comparison with SGR 19354-2154 radio bursts 


Table 1 provides the peak flux density of FRB 20200120E bursts. At a distance of 3.6 Mpc, the 
isotropic radio luminosity of the bursts would be — 10?" erg s^! similar to those of the very bright 
SGR 19354-2154 radio bursts recently detected by CHIME/FRB and STARE2 (The CHIME/FRB 
Collaboration et al. 2020; Bochenek et al. 2020). In 2 yrs of CHIME/FRB observations, we have 
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seen at least three bursts from the FRB 20200120E source. There are other low-DM FRBs within 
the CHIME/FRB sample that are presently under consideration. Careful study of the host galaxy 
candidates in their error regions (which are presently mostly larger than for FRB 20200120E) must 
be done to assert them as extragalactic, given the ever-present possibility of them being in the distant 
MW or MW halo. This work is underway and may be able to constrain the local volumetric FRB 
rate to compare with that inferred from SGR 1935+2154. 

The proximity of FRB 20200120E makes it an attractive target for X-ray and gamma-ray telescopes. 
For a fiducial current high-energy telescope fluence detection sensitivity threshold of 1071? erg cm~?, 
high-energy bursts from nearby sources with energies > 10“! erg s^! should be detectable, and are 
sometimes seen from Galactic magnetars in outburst (see Kaspi & Beloborodov 2017, for a review). 
Moreover, giant magnetar flares with total isotropic luminosity ~ 10/6 erg, like that from SGR 
1806—20 (Palmer et al. 2005), should be easily detectable by the Swift/Burst Alert Telescope (BAT) 
and Fermi/GBM which have flux sensitivity of ~ 10^" erg cm ? s^! in the 15-150 keV band (A. 
Tohuvavohu, private communication) and ~ 10^" erg cm ? s^! in 50-300 keV band (von Kienlin 
et al. 2020), respectively. Unfortunately, Fermi/GBM was either not operational (transiting through 
the South Atlantic Anomaly region) or the FRB location was occulted by the Earth at all but one 
burst epoch. At the time of the 2020 February 6 burst, the FRB was visible to Fermi/ GBM but no 
trigger was reported by the Fermi collaboration. Additionally, the Swift/BAT field-of-view (FOV) did 
not cover the FRB sky-position at the time of FRBs 20200120E and 20201129A. However, the FRB 
location was within the Swift/BAT FOV at the time of FRB 20200718A, and no coincident X-ray 
burst was reported by the Swift collaboration. Therefore, if the FRB source is at 3.6 Mpc, it seems 
unlikely that FRB 20200718A and the 2020 February 6 event were associated with SGR 1806—20-like 
giant flares. 


4.3. Comparison with other repeating FRB hosts 


In contrast with the late-type galaxy hosts of the only three other localized repeating FRBs, FRBs 
121102, 180916 and 190711 (Tendulkar et al. 2017; Marcote et al. 2020; Macquart et al. 2020), M81 
is an early-type spiral galaxy of morphology SA(s)ab (Bosma 1981). Table 3 lists its main physical 
properties. M81 also contains a low-luminosity AGN (Markoff et al. 2008), and is classified as a 
LINER Seyfert (Ho et al. 1996). Heintz et al. (2020) noted that the hosts of apparently non-repeating 
FRBs are typically more massive than those of repeating FRBs. However, M81 would be among the 
most massive FRB hosts known thus far! with stellar mass of 7.2 x 10'? Mo (see Table 3). Lastly, if 
we ignore FRB 190523 for which the host association is not firm (Macquart et al. 2020; Heintz et al. 
2020), FRB 20200120E would show the largest projected offset from the center of its host (~ 20 kpc). 
This would be at odds with the offset distribution of the progenitors of long gamma-ray bursts and 
superluminous supernovae, which are found close to the centers of their respective hosts (Lunnan 
et al. 2015; Blanchard et al. 2016; Heintz et al. 2020; Mannings et al. 2020b). 

If FRB 20200120E is a classical magnetar, it would be surprising to find it at such a large offset 
from the center of its host; known Galactic magnetars are all well within the optical disk of the Milky 
Way (Olausen & Kaspi 2014). This same issue has been noted for other localized FRBs (Heintz et al. 
2020). However, the M81 circumgalactic medium (CGM) is dynamic and rich in gas and metals 
(Chen 2017). Sun et al. (2005) and Smercina et al. (2019) noted the existence of a diffuse stellar 
population embedded in the extended Hi disk where in situ star formation is actively taking place. 
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Table 3. Notable properties of M81, the most promising host of FRB 20200120E. 


Property Value Reference 

SFR (Mo yr!) 0.4—0.8 Gordon et al. (2004) 
Metallicity^ [Z] (dex) 0.03 Kong et al. (2000) 

Stellar mass (Mo) (7.2 +1.7) x 10'° de Blok et al. (2008) 
Effective radius (Reg; kpc) 3.5 Sheth et al. (2010) 

(u-r)o? (mag) 2.7 13(4) Abazajian et al. (2009) 
E(V-B)* 0.26 Kudritzki et al. (2012) 
Absolute r-band mag. (AB) —19.78 -— 

Inclination angle (°) 62 Karachentsev et al. (2013) 
Luminosity distance (Mpc) 3.63 + 0.34 Karachentsev et al. (2013) 
Projected FRB offset from galaxy center (kpe)? 20*3 'This paper 


^ Average metallicity relative to the Sun i.e., log(Z/Zo). However, at the FRB location, the metallicity is 
found to be sub-solar, i.e., [Z] « 0 (Williams et al. 2009). 

b Milky Way extinction is corrected using the reddening map by Schlegel et al. (1998). 

^ Average value of the color excess; at the FRB location, it is likely to be < 0.1 (Kudritzki et al. 2012). 

4 9096 confidence interval. 


Lastly, Frederiks et al. (2007) and Hurley et al. (2010) have argued for the existence of a neutron 
star population, including a possible magnetar, in the CGM of M81. 


5. CONCLUSIONS 


We have reported on the detection of the repeating fast radio burst source FRB 20200120E discov- 
ered with CHIME/FRB. This source has very low DM, 87.82 pc cm ?, though greater than what is 
expected from models of the Milky Way ISM along its line-of-sight. Due to large uncertainties in the 
Milky Way halo DM contribution, it is possible that the FRB source is within our halo. However, 
we find no cataloged Milky Way halo satellite galaxy or globular cluster within or near to the FRB 
localization region that can host the FRB source. Moreover, the presence of a solitary neutron star 
capable of producing FRB-like radio emission at a distance ~ 50 — 200 kpc seems unlikely. On 
the other hand, we identify M81, a nearby grand-design spiral galaxy at a distance of ~ 3.6 Mpc, 
with an angular offset ~ 19’ and chance coincidence probability < 10~?, making it a promising host 
candidate. 

We have shown that the observed extragalactic DM component of the FRB is significantly lower 
than the model-predicted DM contribution from the M81 halo as a foreground galaxy. This suggests 
that the FRB host galaxy is unlikely to be located beyond M81, though the FRB may exist within the 
extended disk of M81. Therefore, if extragalactic, FRB 20200120E is most likely associated with M81. 
M81 is different from the hosts of other known repeating FRBs in spatial offset, stellar population age, 
and local environment. This supports an interesting diversity in repeater host properties that addi- 
tional localizations will help understand. We also found that the FRB localization region contains the 
extended M81 HI-disk and a number of interesting M81 sources including an H 1 region ([IPWK2012] 
31), an X-ray binary ([SPZ2011] 8) and a VLASS source, VLASSIQLCIR J095756.104-684833.3. At 
a distance of 3.6 Mpc, it should be possible to detect prompt multi-wavelength counterparts of FRB 
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20200120E predicted by several FRB models, including the magnetar model. Some FRB models 
anticipate even greater luminosities in high energy bands than in the radio band (Yi et al. 2014; 
Burke-Spolaor 2018; Chen et al. 2020). For example, in the synchrotron maser model of Metzger 
et al. (2019), shock-heated electrons gyrate to produce synchrotron radiation that sweeps through 
the y-ray and X-ray bands, and in some cases, even extends to the optical band on sub-second 
timescales. Additionally, if radio bursts from FRB 20200120E are accompanied by X-ray bursts, as 
was seen in SGR 1935+2154 (The CHIME/FRB Collaboration et al. 2020; Bochenek et al. 2020; 
Mereghetti et al. 2020; Li et al. 2021; Ridnaia et al. 2020), detecting a coincident X-ray counterpart 
seems feasible, and would be a strong test of the magnetar origin of FRBs. Therefore, we encourage 
multi-wavelength follow-up of FRB 20200120E. 
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